PRIME DIVISORS OF THE LAGARIAS SEQUENCE 



PlETER MOREE AND PETER StEVENHAGEN 

Abstract. We solve a 1985 challenge problem posed by Lagarias [5] by determining, 
under GRH, the density of the set of prime numbers that occur as divisor of some 
term of the sequence {x n }'^ } =1 defined by the linear recurrence a; n +i = x n + x n — i 
and the initial values xo = 3 and x\ = 1. This is the first example of a 'non-torsion' 
second order recurrent sequence with irreducible recurrence relation for which we can 
determine the associated density of prime divisors. 



1. Introduction 

In 1985, Lagarias [5] showed that the set of prime numbers that divide some Lucas 
number has a natural density 2/3 inside the set of all prime numbers. Here the Lucas 
numbers are the terms of the second order recurrent sequence {x n }^ =1 defined by 
the linear recurrence x n +i = x n + x n -i and the initial values xq = 2 and x\ = 1. 
Lagarias' s method is a quadratic analogue of the approach used by Hasse [2, 3] 
in determining, for a given non-zero integer a, the density of the set of the prime 
divisors of the numbers of the form a n + 1. Note that the sequence {a n + 1}^ =1 also 
satisfies a second order recurrence. 

Hasse and Lagarias apply the Chebotarev density theorem to a suitable tower of 
Kummer fields. Their method of 'Chebotarev partitioning' can be adapted to deal 
with the class of second order recurrent sequences that are now known as 'torsion 
sequences' [1, 8, 11]. For second order recurrent integer sequences that do not enjoy 
the rather special condition of being 'torsion', it can no longer be applied. In the 
case of the the Lucas numbers, changing the initial values into xq = 3 and x± = 1 
(while leaving the recurrence x n+ i = x n + x n -i unchanged) leads to a sequence 
for which Lagarias remarks that his method fails, and he wonders whether some 
modification of it can be made to work. 

We will explain how non-torsion sequences lead to a question that is reminiscent 
of the Artin primitive root conjecture. In particular, we will see that for a non- 
torsion sequence, there is no number field F (of finite degree) with the property that 
all primes having a given splitting behavior in F divide some term of the sequence. 
It follows that Chebotarev partitioning can not be applied directly. However, it 
is possible to combine the technique of Chebotarev partitioning with the analytic 
techniques employed by Hooley [4] in his proof (under assumption of the generalized 
Riemann hypothesis) of Artin's primitive root conjecture. In the case of the modified 
Lucas sequence proposed by Lagarias, we give a full analysis of the situation and 
prove the following theorem. 
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Theorem. Let {x n }^L be the integer sequence defined by 1 and the 

linear recurrence x n +i = x n + x n -\. If the generalized Riemann hypothesis holds, 
then the set of prime numbers that divide some term of this sequence has a natural 
density. It equals 

1573727 t— r / p 



II ( X ~ ^~t) ~ -577470679956 . 



1569610 

p prime 

Numerically, one finds that 45198 out of the first 78498 primes below 10 6 divide 
the sequence: a fraction close to .5758. 

2. Second order recurrences 

Let X = {x n }^ =0 be a second order recurrent sequence. It is our aim to determine, 
whenever it exists, the density (inside the set of all primes) of the set of prime 
numbers p that divide some term of X. 

We let x n+ 2 = aix n +i + aox n be the recurrence satisfied by X, and denote by 
f = T 2 — aiT — aQ G Z[T] the corresponding characteristic polynomial. We factor / 
over an algebraic closure of Q as / = (T — a)(T — a). 

In order to avoid trivialities, we will assume that X does not satisfy a first order 
recurrence, so that aa = ao does not vanish. The root quotient r = r(f) of the 
recurrence, which is only determined up to inversion, is then defined as r = a/a. It 
is either a rational number or a quadratic irrationality of norm 1. In the separable 
case r^lwe have 

n , ~~n .,, x l ~ ® x , ~ Xi- ax 

x n = ca + ca with c = ^— and c = — . 

a — a a — a 

As our sequence is by assumption not of order smaller than 2, we have cc ^ 0. 
Denote by 

Xi CXXq . . . + 

q = — = — c/c G Q(a) 

xi — axo 

the initial quotient q = q(X) of X. Just as the root quotient, this is a number 
determined up to inversion that is either rational or quadratic of norm 1. The 
elementary but fundamental observation for second order recurrences is that for 
almost all primes p, we have the fundamental equivalence 

p divides x n ^ -c/c = (a/a) n G O/pO q = r n G {O/pOf. 

Here O is the ring of integers in the field generated by the roots of /. This is 
the ring Z if / has rational roots, and the ring of integers of the quadratic field 
Q[X]/(f) = Q(\/af — 4ao) otherwise. The equivalence above does not make sense 
for the finitely many primes p for which either r or q is not invertible modulo p, 
but this is irrelevant for density purposes. 

In the degenerate case where the root quotient r is a root of unity, it is easily 
seen that the set of primes dividing some term of X is either finite or cofinite in the 
set of all primes. We will further exclude this case, which includes the inseparable 
case r = 1, for which q is not defined. 

As we are essentially interested in the set of of primes p for which q is in the 
subgroup generated by r in the finite group (O/pO)* of invertible residue classes 
modulo p, we can formulate the problem we are trying to solve without any reference 
to recurrent sequences. Depending on whether the root quotient r is rational or 
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Problem 1. Given two non-zero rational numbers q and r 7^ ±1, compute, when- 
ever it exists, the density of the set of primes p for which we have 



Problem 2. Let r be a quadratic irrationality of norm 1 and O the ring of integers 
°/Q( r )- Given an element q G Q(r) of norm 1, compute, whenever it exists, the 
density of the set of rational primes p for which we have 



The instances of the two problems above where (q mod r) is a torsion element in the 
group Q(r)*/ (r) are referred to as torsion cases of the problem, and the sequences 
that give rise to them are known as torsion sequences. The sequences {a n + 1}^Lq 
studied by Hasse, the Lucas sequence {e n + £~ n }^ = o with e = 1+ 2 V ^ treated by 
Lagarias and the Lucas-type sequences in [8] are torsion; in fact, they all have 
q = — 1. The main theorem for torsion sequences, for which we refer to [11], is the 
following. 

Theorem. Let X be a second order torsion sequence. Then the set tvx of prime 
divisors of X has a positive rational density. 



Problem 1 in the previous section is reminiscent of Artin's famous question on 
primitive roots: given a non-zero rational number r 7^ ±1, for how many primes p 
does r generate the group F* of units modulo pi (One naturally excludes the finitely 
many primes p dividing the numerator or denominator of r from consideration.) 
Artin's conjectural answer to this question is based on the observation that the 
index [F* : (r)] is divisible by j if and only if p splits completely in the splitting 
field Fj = Q(Cj, r 1 /- 7 ) of the polynomial X^ — r over Q. Thus, r is a primitive root 
modulo p if and only p does not split completely in any of the fields Fj with j > 1. 
For fixed j, the set Sj of primes that do split completely in Fj has natural density 
l/[Fj : Q] by the Chebotarev density theorem. Applying an inclusion-exclusion 
argument to the sets Sj : one expects the set S = Si\ Uj>iSj of primes for which 
r is a primitive root to have natural density 



Note that the right hand side of (3.1) converges for all r G Q* \ {±1} as [Fj : Q] is 
a divisor of <p(j) ■ j with cofactor bounded by a constant depending only on r. 

A 'multiplicative version' of the 'additive formula' (3.1) for 5(r) is obtained if 
one starts from the observation that r G Q* \ {±1} is a primitive root if and only 
if p does not split completely in any field F# with I prime. The fields Fg are of 
degree £(£ — 1) for almost all primes £, and using the fact that they are almost 
'independent', one can successively eliminate the primes that split completely in 
some Fg to arrive at a heuristic density 



q mod p G (r mod p) C F* 



q modp G (r mod p) C (O/pO)*. 



3. Non-torsion sequences. 




oc 



(3.1) 



(3.2) 



5(r)=c r - (1 



1 



)=Cr- C 1 



\F, : Ol 



£(£-1) 
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The correction factor c r for the 'dependency' between the fields Fg is equal to 1 if 
the family of fields {Fiji is linearly disjoint over Q, i.e., if each field Fg is linearly 
disjoint over Q from the compositum of the fields Fg with £ ^ £ . If r is not a 
perfect power in Q*, we have c r = c r . 

It turns out that the only possible obstruction to the linear disjointness of the 
fields Fg occurs when F 2 = Q(\A") is quadratic of odd discriminant. In this case, 
F 2 is contained in the compositum of the fields Fg with £ dividing its discriminant. 
The value of c r is a rational number, and one can derive a closed formula for it as 
in [4, p. 220]. 

For example, taking r = 5, one has F 2 C F5 and the superfluous 'Euler factor' 
1 - [F 5 : Q]" 1 = |§ at £ = 5 in the product U e (l - [Fg : Q]" 1 ) is 'removed' by the 
correction factor C5 = C5 = || . 

It is non-trivial to make the heuristics above into a proof. As Hooley [4] showed, 
it can be done if one is willing to assume estimates for the remainder term in the 
prime number theorem for the fields Fj that are currently only known to hold under 
assumption of the generalized Riemann hypothesis. One should realize that only 
when we consider finitely many £ (or j) at a time, the Chebotarev density theorem 
gives us the densities we want. After taking a 'limit' over all £, we only know that 
the right hand side of (3.1) or (3.2) is an upper density for the set of primes p for 
which r is a primitive root. We have however no guarantee that we are left with 
a non-empty set of such p. Put somewhat differently, we can not obtain primes p 
for which (r mod p) is a primitive root by imposing a splitting condition on p in a 
number field F of finite degree; clearly, there is always some field Fg that is linearly 
disjoint from F, and no splitting condition in F will yield the 'correct' splitting 
behavior in .Fg. A similar phenomenon occurs in the analysis of non-torsion cases of 
the Problems 1 and 2. This is exactly what makes non-torsion sequences so much 
harder to analyze than torsion sequences. 

If (q mod r) is not a torsion element in Q*/ (r), then Problem 1 can be treated 
by a generalization of the arguments used by Artin. For each integer i > 1, one 
considers the set of primes p (not dividing the numerator or denominator of either 
q or r) for which the index [F* : (r)] is equal to i and the index [F* : (q)] is divisible 
by i. These are the primes that split completely in the field F it i = Q(Ci, r 1 /*, q 1 ^ 1 ), 
but not in any of the fields F it j = Q(Cij, r 1 / 1 - 7 , q 1 ^) withj > 1. As before, inclusion- 
exclusion yields a conjectural value for the density 5i(r, q) of this set of primes, and 
summing over i we get 



as a conjectural value for the density in Problem 1. Note that 5i(r, q) is nothing 
but the primitive root density 5(r) from (3.1). 

The condition that (q mod r) is not a torsion element in Q*/( r ) means that q 
and r are multiplicatively independent in Q*. In this case [Fij : Q] is a divisor of 
i 2 j ■ <p(ij) with cofactor bounded by a constant depending only on q and r. Thus 
the double sum in (3.3) converges, and under GRH one can show [9, 10] that its 
value is indeed the density one is asked to determine in Problem 1. 

As in Artin's case, one can obtain a multiplicative version of (3.3) by a 'prime- 




00 



oo 



(3.3) 
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in F* means that for all primes £, we have an inclusion 

(3.4) (q mod p) e C (r mod p) t 

of the ^-primary parts of these subgroups. If we fix both I and the number k = 
orde(p — 1) of factors £ in the order of F*, this condition can be rephrased in terms 
of the splitting behavior of p in the number field 

(3.5) n<*> = Q(^ +1 y/'V/<*). 

More precisely, we have ordi(p — 1) = k if and only if p splits completely in Q(C^fc) 
but not in Q(Qk+i); of the primes p that meet this condition, we want those p 
for which the order of the Frobenius elements over p in Gal(Q((^fc, q 1 ' 1 )/Q(Ct k )) 
divides the order of the Frobenius elements over p in Gal(Q((^fc, r 1 ^ )/Q(Q k ))- By 
the Chebotarev density theorem, one finds that the set of primes p with ord^(p— 1) = 
k, which has density £~ k for k > 1, is a union of two sets that each have a density: 
the set of primes p for which the inclusion (3.4) holds and the set of p for which 
it does not. This 'Chebotarev partitioning' allows us to compute, for each £, the 
density of the primes p for which we have the inclusion (3.4): summing over k in 
the previous argument yields a lower density, and this is the required density as we 
can apply the same argument to the complementary set of primes. 

For all but finitely many £, the fields Q((^fc+i , g 1 ^ ) and Q((^fc+i , r 1 ^ ) are 
linearly disjoint extensions of Q(Qk+i) with Galois group Z/£ k Z for all k > 0. In 
this case the set of primes p with ordg(p — 1) = k violating (3.4) has density 

k 

£- k J2£~\£~ {t ~ 1) '£~ l ) = (£~ k -£~ 3k )/(£ + i). 

Summing over /c, we find that (3.4) does not hold for a set of primes of density 
£/(£ 3 - 1). As the fields 

fl/ = Ufif ) Q(0-,r 1 /'V /< ") 

for prime values of £ form a linearly disjoint family if we exclude finitely many 'bad' 
primes £, the multiplicative analogue of (3.3) reads 

(3-6) S(r,q)c q , r - J] U - . 

I prime 

As is shown in [9], the 'correction factor' ^ f IS rational number that admits a 
somewhat involved description in terms of q and r. In practice, one finds its value 
most easily by starting from the additive formula (3.3). 

In the situation of Problem 2, the arguments just given can be taken over without 
substantial changes from the rational case when one restricts to those rational 
primes p that split completely in O. Writing K = Q(r) and 
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we find that the density of the rational primes p that are split in O and for which 
we have q mod p G (r mod p) C (O/pO)* equals 

OO OO / .x 

i=l j=i ™ • 

as in the case of Problem 1, one needs to assume the validity of the generalized 
Riemann hypothesis for this result. By (3.7), the computation of 5 sp n t amounts 
to a degree computation for the family of fields {Fij}ij. In fact, because of the 
numerator fi(j) in (3.7) one may restrict to the case where j is squarefree. As in 
the case of Problem 1, one finds (under GRH) that the split density equals 

(3-8) * B piit(r,«)c+ p . II ( 1_ F^t) 

£ prime 

for some rational number c+ r . The next section provides a typical example of such 
a computation. It shows that the value of c+ r is not as simple a fraction as the 
analogous factor c r in (3.2). 

For the rational primes p that are inert in O, the determination of the corre- 
sponding density 5 ine rt is more involved than in the split case. The group (O/pO)* 
in Problem 2 is now cyclic of order p 2 — 1, and (q mod p) and (r mod p) are elements 
of the kernel 



(3.9) k p = ker[N : (O/pO)* — ► F p ] 

of the norm map, which is cyclic of order p + 1. In order to have the inclusion (3.4) 
of subgroups of k p for all primes £, we fix £ and k = ordi(p + 1) > 1 and rephrase 
(3.4) in terms of the splitting behavior of p in the quadratic counterpart 

nP = K(Q k+i y/e k ,q^ k ) 

of (3.5). Let us assume for simplicity that £ is an odd prime, and that K is not the 
quadratic subfield of Q(Q). Then the requirement that p be inert in K and satisfy 
ord^(p+ 1) = k > 1 means that the Frobenius element of p in Gal(K (Qk+i) / Q) is 
non-trivial on K and has order 2£ when restricted to Q((^k+i). Let Bk C K(Qk+i) 
be the fixed field of the subgroup generated by such a Frobenius element. Then Bk 
does not contain K or Q(0), and B}~ C K(( e k) is a quadratic extension. Let 
be the non-trivial automorphism of this extension. Then cr/. acts by inversion on 
(gk, and the norm-l-condition on q and r means that cr*; also acts by inversion on q 
and r. The Galois equivariancy of the Kummer pairing 

Gal(K(C £k ,r 1/£k ,q 1/ik )/K(C £k )) x (q,r) — 

shows that the natural action of at on Gdl(K(Qk, r 1 /^, q 1 /^)/ K(Q k )) is trivial, 

(k) 

so ' is abelian over It is the linearly disjoint compositum of the cyclotomic 
extension Bk C Sfc(O fe + 1 ) an d the abelian extension 



T~> , — T~> / 1/^ , -1/^ 1/^ , -l/£ h \ 
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For almost all £, the group Gal(of } /B k ) is isomorphic to Z/2£Z x (Z/£ k Z) 2 . 

Just like in the rational case, we want those primes p that have splitting field 
Bj~ inside K(Qk+i) and for which the order of the Frobenius elements over p in 
Gal(Sfc(q |1 ^ + q~ x ^ )/Bk) divides the order of the Frobenius elements over p in 
Gal(S/ c (r 1 /^ -|-r -1 /^ )/Bk). By the Chebotarev partition argument, we find again 
that for a 'generic' prime £, a fraction £j{£? — 1) of the primes p that are inert 
in O violates (3.4). Here 'generic' means that £ is odd and that has degree 
2£ 3k (£ — 1) for k > 1. Under GRH, one can again deduce that the inert density 
<5inert((/7 ?~) equals a rational constant c~ r times the infinite Euler product occurring 
in (3.6) and (3.8). 

In general, there are various subtleties that need to be taken care of in the 
analysis above for £ = 2, when 2-power roots of q and r are adjoined to K or 
K((2 n )- We do not go into them in this paper. In the example in the next section, 
we deal with these complications by combining a simple ad hoc argument for a few 
'bad' £ with the standard treatment for the 'good' £. 

4. The Lag arias example 

We now treat the explicit example of the modified Lucas sequence which is the 
subject of the theorem stated in the introduction. The roots of the characteristic 
polynomial X 2 —X — 1 of the recurrence are e = 1+ 2 V ^ and its conjugate e = 1 ~ 2 V ^ . 
The initial values xq = 3 and x\ = 1 yield an initial quotient q = j^ff of the 
sequence. As 7Tn = 1 — 3e G O = Z[e] has norm —11, we find that we have to solve 
Problem 2 for 

7Tn -k\ x e 2 

q = — — = and r = - = — s . 

7Tn —11 e 

We set K = Q(e) = Q(\/5) and Fij = K(Qj, r 1 / 1 - 7 , q 1 ^) as in the previous section. 

4.1. Lemma. Fori.j G Z>i we have [Fij : Q] = 2 1 ~ t i 2 jLp(ij), with 

#{<i G {4, 5, 11} : d\ij} if i is even; 
* = = { e {4? 5} : d\ij} if i is odd and j is even; 

#{d G {5} : <i|r/'} i/z? is odd. 

Proof. As K = Q(v^5) is the quadratic subfield of Q(Cs) 5 the field K(Qj) has 
degree 2ip(ij) over Q if 5 does not divide ij and degree <p(ij) if it does. 

As q = is a quotient of two non-associate prime elements in O and e is 
a fundamental unit in O, the polynomials X 1 — g and X y — r are irreducible in 
K[X] for all i,j6 Z>i by a standard result as [6, Theorem VI. 9.1]. Moreover, the 
extension K C K(q |1 / 1 ) generated by a zero of X* — q is totally ramified at the 
primes of K lying over 11, whereas the extension K C i^r 1 / 1 - 7 ) generated by a zero 
of X lJ — r is unramified above 11. It follows that K C K(q 1 ^, r 1 /*- 7 ) is of degree i 2 j 
for all i, j G Z>i. 

The intersection K{q 1 / 1 , r 1 /*- 7 ) ri-KT(Cij) is contained in the maximal abelian sub- 
field Kq of K{q 1 / 1 , r 1 / 1 - 7 ), which equals 



K{^/q, y/r) = If (v— IT, C 4 ) if z is even; 
ifo = ■{ K(y/r) = K ((4) if z is odd and j is even; 
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One trivially computes K D K(Qj), and the lemma follows. □ 

We will need the preceding lemma only for squarefree j. In this case, we simply 
have t = #{d G {5} : d\ij} for odd i. 

If we substitute the explicit degrees from Lemma 4.1 in (3.7), we find that the 
split density for our example equals 



OO OO , .n 



where Uj is as in Lemma 4.1. If we set 



oo oo 



m\i mn\ij 

then the expression above may be rewritten as 

25 sp iitSl ! i + S2,2 + S2,5 + S2, 11 + £2,10 + £2,22 + £2,55 + 52,110 + (Si, 5 — S2,5)- 

It is elementary to show [9, Theorem 4.2] that S m , n is the rational multiple 

S -p-r -p 4 yi P 3 +P 2 



1 p s — p — 1 p 3 — p — 1 



p\n P\ m 

of the universal constant S = Y[ p prime (l — p^rzi) for m, n £ Z>i. Simple arithmetic 
now yields the value 

712671 

(4 ' 2 > ^"'1569610 26151 

for the density (under GRH) of the primes p = ±1 mod 5 dividing the Lagarias 
sequence. Numerically, one finds that 20416 primes out of the 78498 primes below 
10 6 are split in K and divide our sequence: a fraction close to .2601. 

For the inert primes of K, which satisfy p = ±2 mod 5, we have a closer look at 
the 'bad' primes 2, 5 and 11. As for the rational problem, we define the extensions 

for primes I and note that, by Lemma 4.1, the family consisting of the extensions 
O2O5O11 and {0^)^2,5, 11 of K is linearly independent over K. 

Our first observation is that for the inert primes p, the order p + 1 of the group 
k p in (3.9) is never divisible by 5. Condition (3.4) is therefore automatic for the 
prime t = 5, and we can disregard the splitting behavior of p in ^5. 

We next observe that for inert p, the element r = — e 2 satisfies 

(4.3) r^+ 1 )/ 2 = (-l)^- 1 )/ 2 modp. 

For primes p = 3 mod 4, this shows that (r mod p)2 is the 2-Sylow subgroup of k p , 
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primes congruent to 3 mod 4, which form a set of primes of density 1/4, have the 
correct splitting behavior in the extensions fig for £ ^ 2, 5, we are dealing with 
a linearly disjoint family and find (under GRH) that the set of these primes has 
density 

4 II V P-l) 4 5 119 ' 



We next consider the inert primes p = 1 mod 4. For these p, the congruence (4.3) 
shows that (r mod p) has odd order in k p , so (3.4) is satisfied for £ = 2 if and 
only if the order of q = — n^/ll in k p is also odd. As k is a cyclic group of order 
p + 1 = 2 mod 4, the order of q = (q mod p) is odd if and only if q is a square in k p . 
Let x G O/pO be a square root of q. If a; is in k p , i.e., if a; has norm 1 in F p , then 
its trace x + 1/x is in F p , and we find that 

1x2 . - 1 1 
(xH — ) = 2 + q + - = mod p 

x q —11 

is a square modulo p. If x is not in k p , then x has norm —1 in F p and 

^2 o - 1 32 • 5 
(x J = 2 — g — - = mod p 

x q —11 

is a square modulo p. As 5 is not a square modulo our inert prime p, we deduce 

(4.4) (q mod p) has odd order in n p —11 is a square modulo p. 

If p satisfies the equivalent conditions of (4.4), then p is a square modulo 11 by 
quadratic reciprocity, and we have 11 \ p+ 1. It follows that in this case, (3.4) is 
satisfied for £ = 2,5 and 11. Thus, the set of inert primes p = 1 mod 4 satisfying 
the quadratic condition (4.4) is a set of primes of density 1/8, and the subset of 
those p that have the correct splitting behavior in the extensions Qi for £ 7^ 2, 5, 11 
has (under GRH) density 

I TT (1 I A _ 1 7 124 1330 
8' 11 I ~ £337 )- 4 ' 5 ' Y19 ' 1319 ' 

£/2,5,ll v 7 



Adding the fractions obtained for the inert primes congruent to 3 mod 4 and to 
1 mod 4, we obtain 

61504 

<Wrt = S w 0.3159598798268. 

112115 

Numerically, one finds that 24781 primes out of the 78498 primes below 10 6 are 
inert in K and divide our sequence: a fraction close to .3157. 

The sum <S split + 5 inert is the value ( + tehis) ' ^tSUS ' S mentioned in 
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